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tion of optimal solutions. Revised arguments are given to 
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iN TROB@CT ION 


In any inventory model, one of the main objectives is 
momavotc bOmnmeversteck and Out-of=-stock Situations. The 
latter is especially critical in military supply. In lot 
Size - reorder point models with stochastic demands, a 
decision has to be made as to when to reorder and how much 
the reorder quantity should be. The goal, then, is -to 
achieve in some sense an optimal reorder level r and an 
optimal reorders quantity Q. One way to accomplish this goal 
is to define a cost function K which then acts as an objective 
function to be minimized by choice of Q and r as variables. 

Hadley and Whitin [{Ref. 1] present several lot size - 
reorder point models. The one of interest in this paper is 
called a backorder model wherein, if demands occur when the 
system is out of stock, they will be backordered with some 
PeemeyecGst,—to be called the backorder cost. Customers 
then must wait until all orders are delivered te the inventory 
manager. Associated with delivery is a delay in time called 
feaeetime win this paper we consider only the case where 
lead-time demand is assumed to be normally distributed. 

Mathematically, the backorder model of interest here has 


the objective function, 


K(Q,r) = B+ ref ere + Beer o(F + oF) 





where 

A is the set-up cost (cost of placing an order) 

tT is the backorder cost per unit backordered 

XA is the average annual demand 

I is the inventory carrying charge 

C is the unit cost of the item 

u and o are the mean and standard deviation of the 
normally distributed lead-time demand random 
variable 

6(z) is the standard normal density function 

O(2) is ftsmeompetementary cumulative distribution function, 


CO 


i,en;meecneze = Jf dfx) dx 
Zz 
iSi@peeetemea tea the average annual cost and 1s composed 
of three different types of cost: the ordering cost, the 
cost of carrying inventory and the cost of backorders. These 
costs are represented, respectively, by the first, second and 
Emad term in the eee nomeece 1) dia (1). 


For ease in the algebra, (1) will be rewritten as 


K(Q,z) = B+ rctZ + oz] + MF [o(z) - 20(2)] (2) 
Whete= 2 — — 


Heog@leywendewaltin [Ref. 1, pp. 165 ££.}] point out that 
the model just defined is appropriate when the expected number 
of backorders is negligible. They then propose a solution to 


the optimization problem by claiming that K is convex in Q 





and r (or z) and therefore any solution obtained by setting 
the partial derivatives of K equal to zero will determine a 
global minimum. They provide an iterative scheme for such 

a solution along with the assertion that the solution, when 
Me exists es Unvague. Finally, they discuss the fact that 

a solution will not always exist, pointing out that this 
"anomaly" arises because a backorders term was omitted in the 
expression for the carrying cost. Since the model is to apply 
only when such a term is negligible, they assert that this 
Situation will never occur when high backorder costs are 
involved without being explicit as to what "high" means. 

To verify these results, the reader is directed to solve 
Peso atccnommeyroccises, the first of which is to show that 
PemruneevOn Involving the backorder cost (the third term) 
Deel say 0(0,r), 1S convex in QO and r. It would then 
follow immediately that K is convex in Q andr. But it was 
Memnted Out as early as 1964 by Veinott [Ref. 2] that J is 
not convex. In 1969, Brooks and Lu [Ref. 3] addressed the 
Same problem and derived a general result which, when applied 
Eemthe normal case, shows that J is convex for all rou, 
the mean lead-time demand. 

However, these results still leave the main issue 
unresolved. In the first place, solutions can be obtained 


where x < wet More importantly, the real isSue is whether 


the difference r-u is sometimes referred to as a 
safety level. In that case, a negative safety level simply 
Means 1t) 1s more economical to tolerate a backorder position 
on the average than to protect against that position. 








or not K is convex, and, in any case, how its minimum value 
iowcMeracterizea.. | (he purpose of this thesis is to clarify 
this issue and characterize the solutions to the problem in 
a more explicit manner. 

Bectlogel. wil show that in fact K is not convex and, 
in addition, will exhibit an example in which the first 
partial derivatives a and - vanish at two different 
points. Thus, the solution is not unique as claimed by 
Hadley and Whitin. Section III will deal with the character- 
ization of the solutions and will show that under suitable 
conditions there are always two distinct solutions, one of 
them being a@ minimum and the other a Saddle point. Conditions 
for the existence of a solution will be given in terms of the 
Set-up cost A. Section IV will discuss the application of 
the results in Section III to practical inventory problems 


along with the iterative scheme given by Hadley and Whitin to 


solve for a solution. 





10d BA 


COUNTERE XAMPLES 


As pointed out in the introduction, the backorders term 


imeche cOShefunction is not 
and Whitin [Ref. 1, Problem 
still be the case that K is 


Maat this is not so will be 


convex as asserted by Hadley 


S07). 221}%. “Even so, 1t might 


convex as asserted by the authors. 


demonstrated in part A. On the 


other hand, even though K is not convex, it might still be 


the case that the partial derivatives vanish at a unique 


point and that point determines a minimum. 


IM oGEE Wie)lake Sle javere 


sO will be shown in part B and a characterization of solutions 


will be given in the next section. 


Pewee el VEX TY CONDITION FOR 


K 


The function K 1s convex if and only if its Hessian 


Matrix 2 1S positive semi-definite [Ref. 


4]. Since K is a 


function of two variables, the matrix 2 will be positive 


semi~definite when its diagonal terms and its determinant 


are nonnegative. 


For the case at hand, using (2), the Hessian matrix for K 


1s given by 


> 


2 


.@) 


TAO (Zz) 


(A + amn(z)) 


TAO(Z) 


TA > (Zz) 
oQ 
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where n(z) Ooo(z) - 026(z) 


Or ny, Gudz) (3) 


Oz) = 20(z) . (4) 


where 0 (4) 


in Appendix A it 1s shown that u(z) > 0 for all z 
Hence, for nonnegative set-up costs A, it is clear that the 
Giagonal terms of 2 are positive. Thus K is convex if and 


Peyeet det 2 se) 0 . 


Now, 
2 De 
ee 277A mee i ee TA a2) 
oQ 0 
emcethus det > 0 if and only if 
2 Drage 
4 4 nile 
oQ 0 


Multiplying the above inequality by of > 0 and dividing 


by me > Umea yes 


2942) (y + tn(z)) - 10° (z) cameo 
or 2 o(z) + 216% (z) =m (2)O( 2) = 16° (z) > 0 (5) 


Defining V(z) to be equal to the left-hand side of (5), then 


the following summarizes convexity. 


li 





Meee ne esiunetaon K issconvex if and only if det 2 > 0 


or equivalently if and only if WN ¢e |) ee 18 


REMARK: This result holds for any Q > 0 and it is shown 


in Appendix A that V(z) > 0 whenever z > 0 


A counterexample to the convexity of K is then attainable 
from an example provided by Hadley and Whitin [Ref. l, p. 173]. 
There the parameters chosen are 7 = 2000, A = 1600 , C#=Q 50, 
I= .20 , A = 4000 , py = 750 , oo = 50 . For such a case 


V(z) becomes 


Tee et 000e2 (2) = 4000z26(z)o(z) - 200007 (z) 


and for 2z2= -l (equivalently, r = 700) 


V(-1) Beam ce 2o4.256 + ° 894 >3784 - 1415.57138 


-~328.21698 < 0 


Thus, K cannot be convex. 


Bi. COUNTEREXAMPLE TO UNIQUENESS OF THE SOLUTION 
DEFINITION. (Q,z) will be called a solution whenever Q > 0 


and z satisfy the following equations: 


AA, FC _ WAGulz) 


a = — + = = — = ee Ce Ct (6) 
dQ he 2 oki he Z ok 
_ 0K _ _ 71A4O(z) 
0 = a GiC ca )o (7) 


2 





It is convenient to write (6) and (7) as: 
2X(A + Tn(z)) = 2A(A + Toul(z)) (8) 
C= Ic fC 


OZ) = Th (9) 


and 


COUNTEREXAMPLE. Let 7 =A =o021, I= 4.1, C= 2.42 
= _ il 2 a S 

and u= 3. Let A= Url) (2242 Cat. 05) OOS) ee hess 

mmere u(z) is defined in (4). Then there are two distinct 


solutions 
%(z,) 
is,) (Q,,2,) = ( 1 ew 24) where 24) = .05 


$(z.) 
— Z c= = 
2) (0,125) = ( (SL) (2242) ' Z>) where 1G Les Zo < -1.80 


To see that (Q, 724) St a Gonntton, 


67 (.05) 
ak _ 2 ena 7 809)! ay (2.42) |) ) (v.05) 
80,2, 67 (05) : 7 (.05) 
(.1)7(2.42)° Gy ea 


_ 0) (2.42) , 6.107 (2.42) 20.05), (1) (2-42) _ (192 (2.42) 20.05) 
: 6 (05) : 6 (.05) 


0 . Thus, equation (6) is Satisfied. 
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aK _ _ (1) (1)6(.05) _ _ 7 
_ , pe eee (6 -1)1(2.42)) = (.1) (2.42) = 0 - 
ie] (ene. 42) 


Hence equation (7) is satisfied. Therefore, (Q) +24) is a 


solution. 
Now, let S(z) = : ia - u(z) Since S is 
: pam) (2.42) : 
clearly a continuous function and 


S(-1.81) < .0997 < A < .1060 < S(-1.80) 


there must exist some z. Such that -1.81 < Zo < -1.80 


2 
and S (25) =A. The verification for (OQ, 425) 1s done ina 
Similar manner. 
1 [ 6° (2) Bony) 

OK ee Pee 322) 2 fi 2.42) (1) (1) (1) (22) 
Bo,» 2), o* (22) g o° (22) 

2 2 Z 2 

(or (2.22) (ab) (2.42) 


Meno) (2,42) (22) | (12.42) _ 1) 2 (2.42) ua) 
o* (22) (22) 


=0 . Thus, equation (6) is satisfied. 


OK _ eC ian 22) 7 7 
al _* Calinig42) = Se (.1) (2.42) - (.1) (2.42) = 0 
2'= (.1) (2.42) 


and hence equation (7) is satisfied. Therefore, (Q,+2,) 1s 


also a solution. 
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To illustrate these results, a graph similar to the one 
given on p. 171 of Ref. 1 is reproduced in Fig. 2-1 for this 
smecieal case. In this graph r= oz+u1 is expressed as a 
function of Q for each of Equations (8) and (9). It shows 


that the two curves intersect at more than one point. 


Eq(9) | \Eq(8) 


(Q; Z| ) 


m_—w =~ em ans ww 8 = ot Ge —_—=_ ene == ws af ew es 2p ete eee oD 


te 
ne) 
O 
> 
O1 


rg ea 2= 1% 


les 





i Pe 
The value A = IT) (242 6°(.05) - u(.05) was by no 


means chosen by accident. AS it will be seen in the next 
section, the set-up cost A plays an important role in the 
eigtpacterlzatilon Of tne Solutions, and it will be shown that, 
under appropriate conditions of the other parameters, there 


are always two distinct solutions, one of which is a minimun. 
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TILT. CHARACTERIZATION OF MINIMUM SOLUTION 


Let a set of parameters be given. To simplify matters, 


define 


mh 


5TC aA oh - TO0uU(Zz) (150) 


S(z) = 


Peeseries Of lemmas are convenient to derive the main result 


in a theorem. 


LEMMA 1. Let (0*,z*) be a solution. Then, 


(a) A = S(z*) 


(b) V(z*) > 0 if and only if ¢(z*) > <2 and v(z*) = 0 
mance only 2f o(z*) = =o ; 


Peeotma). By definition of a solution, Q* and z* must 


Satisfy (8) and (9) of Section III. Thus, 





Came mou, Z*)) =. (Q*) 2 2 (1) 297 (24) 

— <C,; — [2D 
(IC) 

One 
(nA) 7 2 
2\A + 27hou(z*) = —zE- O (2%) 
pan = LONE 92 (oe) - amdoulzt) 
= <e do (Zz TAOU(Z 
2 
A = Sie ometz) = Tou(z*) = S(z*) 


tLe 





PreoOot (b). 





V(z*) = ES 9 (z*) meecitowz) (o(z*)—-z*o(Z*)) - One) 
= of Ce ere ro(z*)u(z*) - 10° (z*) 
= 2$(z*) (2 + Tu(z*)) - nO* (z*) 
TX 2 2 
= 2$(z*) (sag? Zena )tmu(zZ*)) — wo (z2*) £Erom part (a). 
= TA g (ney 62 (ae) - 162 (2%) 
= Co Zz zZ Z 
= Td 32 (o% x) = LEG 
= TGs (Za lotz ae! 
mn 2 
Since iGo 6° (z*) > 0, the result stated in (b) follows 
immediately. 


OF. LD < 
It is not always true that there exists a solution 
for any given set of parameters. The following lemma discusses 


conditions for which K cannot be minimized. 


tee 2 eeeGiven A > 0 and 7, A, I, C, 6 all positive, a 


necessary condition for the function K to be minimized is 
reo 
that Th < O00). 
Proof. The proof is divided into two cases. 
ICo 





Case (i): Suppose 1 movomanan. 15 Minimizea at (O*,z*) . 
iienmene Partial derivatives muSt vanish at (Q*,z*). But 
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> 
V(z*) = tec 0” (2*) [6(z*) = =o) < 0 





mec OZ) gee and sO {0*,z*) 1S a saddle point. This 
mor a IcOntradlculon, 


Case (ii): Suppose =e = 6(0) and K iS minimized at (Q*,2*). 


Again, the partial derivatives must vanish at (Q*,z*). Now 
if z* # 0 then V(z*) < 0 and hence (Q*,z*) is a saddle point. 
Thus it must be the caSe that z*=0. But then, by Lemma l, 


part (a), 


2 
UE ie 
A = 5TC ® (0) Tou (0) 


TA 
= SIC = TO 6 (0) 


io} 1 
a: eG OK) 
weno. 
me ot) 7°98" 


= TO ( 


ana 84" (0) | 
d (0) 


Since 1 —- 84-(0) < 0, A < 0 contradicting the fact that 


dee Ue 
BS on Ginn a sR Oxe) 
Thus K cannot be minimized unless ae on GoD ae 
OnE. D. 
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LEMMA 3. Suppose aco —O.0) bet Z > 0 and aa < 0 


O 
. | = Be 1 le) 
be the unique points Such that  9( Zo) = o(z,) a 





Then, with S(z) defined in (10), 
(a) S has a minimum at Zz where Be - S(z,) <0 
(b) S has a maximum at hs where es = aco) > 0 


(c) S is strictly concave increaSing on the interval [oir ae) 


ay) 
Proof. S(z) = Le ae - 710u(z) 
m2 
S'(z) = - to 0 ¢2) ez) = CER A) 
1 
ae aa o(z)@(z) - T0(-6(z)) (see Appendix A) 
2 
= 7106(z) - = o(z) (2) 
nd ICo 


= ato? 2) a OZ) 


Since ume @é(z) > 0, the first derivative S'(z) vanishes 


Omlly at =i and Ae Also, S'(z) < 0 for -Z, <zZ < Zo and 


cz) Oe tor Zz < ao 22 2a that is, S is decreasing 
in the interval (=2e2,) and increaSing in the intervals 


eo, —Z.) and (257%) . 
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Z | 1m) 


s"(2) = - TA 6(z) AS - o(z)) + BA ote) (2002) 
Caen 2 2 
_ tT AO~ (z) Ri Azgi2 OZ) 
> COA) a 


2 
_ Tt Ad(z) = 
= Sa a + z0(z)) TOd (Zz) 


T0o(z) [mee ole et = 20(Z)) - 1) 


ae 
Thus, at Zos 


S"(25) = 0$(z,) [gee (ASS + 206(2,)) = 1 


TA 





= 79$(z,) [1 + yar 25%(2z,) - 11 
2 
= EA 256(z,)0(2z,) > 0 (11) 


Therefore, S haS a minimum at Zo: 


To show that Ae = S(z.) < 0, Suppose, to the contrary, 


eaat S{(Z,) > (ene etnce S'(z) > 0 for z > Zot 


S(z) > S (2) Expo eon eall 2g > 7A and hence lim S(z) > QO. 





Z-ro 
2 : ae 2 
Bac ham Sz) = lim sTC 6"(z) - To limu(z) = OO (see 
Zrow Zero 2>0 
Appendix A). Thus, Ne = S(z)) < 0. This together with 


(11) proves part (a). 


Zl 





Now at -Z_, 


es, ee, eg 
tea ‘TA ~ Zo%'52)) - 


ee) Too (-Zz.) [ 


TA 
Eo 


To ¢(-z,) [1 - Zoic) - 1] 


md 
ne 


= 72 UE) rainy Oras (12) 
Therefore, S has a maximum at stan 
To show that AG = S(-z)) > 0 , again suppose that 
MZ.) “SUeemotmce tim s(zZ) = 0 , there exists some z' 
cimem that S(z') > s(-2,) , for otherwise S(z) < SS Ze) OLE 
Beez amples that lim S{(z) < S(-Z)) < 0. But the existence 
Semsuch zz ' MeWeeecicts the Pace peat ols Maximized at TZ. 
Tis, A, = S(-z,) Eames resule coupled with (12) proves 
Part (b). 


Part (c) can now be easily proved by noting that for 


ZS an <5 Os. 


o(z) + 2(z) < $(z) < o(-2,) = => 


hence ae Woe 2o(z)) < 1 
or mes jonz) + Zo{z)) - 1 < 0 


and therefore S"(z) = nog (2) [pe Cee ( 7) ) — 1io< 0 


Thus, S is strictly concave increasing on the interval (S24) 


The proof of Lemma 3 1S now complete. 
OPM ey oF 
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REMARK: The minimum and maximum in Lemma 2 are clearly unique. 
fe a cen a | 


Moreover, lim S(z) = -2e (see Appendix A). A graph of S(z) 


Zraa 


corresponding to the example of Section II is given in Fig. 3-l. 


S 
er, oes ee Sa anses A. 
Kw 
> 
Lo- == A! ee enna 1 Z 
Fig. 3-1. 


The following lemma shows the importance of A in the 
characterization of a solution. 


iL O16) 


LEMMA 4. Suppose = oO) econ 0 1S given. Then 





ie 2 eeomasolution if and only if A = S(z2*). 


Proort. The necessary part has been proved by Lemma 1, part 


(a). The sufficient condition will now be proved. Suppose 
2 
= 5 ee & ee eA 
ie eZ, ) seq 8 (2*) Tmou(z*) . het Q we 


Then, it can be easily shown by Substitution of A and Q* 
into equations (6) and (7) that (Q*,z*) is a solution. 


ORE DY 
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REMARK: It should be observed immediately from Lemmas 4 and 3 
that a necessary condition for a solution to exist is that 


0) Sos Bee The following theorem summarizes all the results. 


LEG 
TA 


and A, ? oe, ~etenwtoe §0a< A < Ag there exist exactly two 





GAEOREM. Giwen tT, A, 1, C, o all positive and < o(0) 


solutions (Q)*,2,*) and (Q5*,Z5*) Suen. tbat Z4* > Zo* ; 


(Q,*,2,%) is a minimum and (Q,*,Z5*) 1s a saddle point for kK. 


Proof. Since S(z) is strictly concave increaSing in the 


interval (—Sie2 ale and strictly decreasing in the interval 


Me yt. 2 5) with hl S(z) =-©- , S(-2,) = Aa ZOO sand 
S(z,) = ee <0, the function S(z) must have two roots 24 
and Zo where Z, < 725 and “Zn < 25 < 2 Eis LOrmany, 


mesouch that 0 < A< Ags there are two distinct values 


* - 
Zy and Zo * such that Z1 < Za * 9 8 RS Z4* < Zo and 


TAd(Z,*) 
A= S (z,*) = S(Z5*). By merely defining Q,* = Se ae 
TAO (Z5*) 

a ee ; * ¥ 
and Q.* iG , it follows by Lemma 4 that (O4 1Zy ) 
and (Q0,*,Z5*) are solutions. But, oe = z4* Saez implies 
that o(z,*) > <=, and hence by Lemma 1, part (b), V(z,*) SO». 
Thus, (0)%*,2,%) MnemazZes Ky Similarly, Zo* S48 implies 
that $(25*) < a8 and hence by the same lemma V(z5*) < 0 


Thus, (Q5%*,25*) is a saddle point. 


0O.E.D. 
REMARK: Bor A= A, = ace , the only solution possible is 
determined by 2z* = -Zy - But telem,~ -adiZ *) = aoe and thus 
V(z*) = 0, so that no conclusion can be drawn. 
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PV Coreen wt TONAL ASPECTS 


Before attempting to solve a given problem involving a 
set of parameters, some simple tests should be performed to 
see whether or not there is any solution. 


1) The first test is ee < 6(0). 





2) If this test is passed, then find mae and A, = S(-z.). 


Two cases could then occur. 
Case (i): SS a Cee Only solution then is -Z, and Wate) 0 
and it has been pointed out that no conclusion may then be 


drawn. 


Case (ii): A, geome te, OS A < Bea then there exist exactly 
; * * * * 7 * ~ rd 

mwo Solutions (Q, 12 and (Q. 1 Zo Pewien Zo* < Zo < 2,* - 

As shown, (Q,*,2,*) is a minimum and (Q5%*,25*) is a saddle 


then again the only Solution is -z. and 


Meant. If A=A 5 


oO! 
no conclusion can be drawn aS in case (i). 

Hadley and Whitin give an iterative scheme for solving 
oie (Oe yan) Ref, 1. This scheme seems to converge always 
to the minimum value and never to the Saddle point. But there 
is a Simple explanation for this. The scheme begins by choosing 
ietially a Wilson Qn = ips for the value of Q to be 
used in Eq. (9) (See Fig. 2-1) and QW. < Q,*. The iteration 


then proceeds to choose a value of r corresponding to OW from 


Eq. (8). With that value of r a new value of Q is selected 


ZS 





from Eq. (9) and it will always be the case that Q < Q,* 
in that selection. Proceeding in this manner, the pairs 
(Q,x) converge to the solution (Q,*,r,*) and cannot yield 
(Q5*,r5*). Thus, when a minimum exists, this computational 
scheme may be relied upon. 

However, in any case of doubt, a Simple test is furnished 
by checking whether or not 6(z*) > 2 fOr ther SOLUtLON - 

It has been pointed out that the boundary case A = IN 
causes problems in that no conclusion may then be drawn. This 
presents no difficulty for practical uses, however. After 
all, an optimal solution may be given for any A arbitrarily 
close to A. : In an application it would be a rare circumstance 


that produced a set of parameters in which A (typically 


measured in dollars and cents) would be exactly Ao: 
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Vee=CONCLUS ION 


While the results of this thesis only apply to the case 
Of normal lead-time demand, they are of more than mere 
academic interest. In fact, that is a common assumption made 
by Hadley and Whitin along with users of the formulas and 
schemes given by them. And it is a happy coincidence that 
the algorithms given are indeed valid even though the 
assumptions under which they were derived were replete with 
errors. 

Nevertheless, it iS a recommendations of this thesis that 
further study into other types of lead-time distributions be 
made. One of the key features of the reSults derived here 
was the examination of the Hessian matrix with the observation 
that its positivity did not depend on Q. This will still be 
true with a density other than normal and so much of the same 
general technique ought to apply. 

One of the features of this thesis has been to amplify 
and clarify completely the minimum solution. Conditions are 
given for existence,and formulas are given along with simple test 
criteria when existence is assured. A by-product of the 
results is the fact that negative safety levels do indeed 
occur as the unique Optimal solution. In fact, a simple 
examination of a graph like S(z) will show immediately and 
Samplify those values of the set-up cost A for which such 


solutions are valid. Since the solution in Such cases is 





unique, certainly negative safety levels should not be ignored 
or tampered with (truncation toy for example) in an application. 
It should also be noted that the solutions are characterized 
independently of u. 

AS a tinal recommendation, this thesis was restricted to 
the so-called backorders model where a customer must wait for 
delivery when the inventory position is out of stock. There 
are, however, other models such as the "lost sale case" 
discussed by Hadley and Whitin to apply in this situation. 

An analysis Similar to that applied here is Surely called for 
so that resultant solutions are validated and characterized. 
For the authors continually appeal to convexity for drawing 


conclusions in those cases. 


ZG 





APEENDEX A. 


As in Section II, 


VARIATION OF V(z) 


ZrO 


2h (a + mh(z)) HAet2) 
Q fe 
Q = 
O0Ad(z) TA d (z) 
Q ie 
TAd(z) . : ae 
The diagonal term a7. 6 Gal Powmecetuaimly POSITIVE. 
mizye= 0(d(z)) = zo(z)) = ou(z) as defined in (4). 
u'(z) = -z26(z) - (@(z) - z¢(z)) 
= -6(z) < 0 Tape wre lel) we, 
Mz O12) > 0 Form val. Z 
Cleamliy, slim u(z) = lim (¢6(z) = z6(z)) = 0 + © = +a 


bm a= OO 





and 
lim u(z) = lim ($(z) - z6(z)) = - lim z4(z) 
7->0o 7->00 7-00 
= - lim 
ae Oz) 
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Applying L'Hospital's Rule repeatedly, 








2 
et eno eee © fz) Cg. = 26(2z) (2) 
rome Te Pe 
Ra oe, 2 
$~ (z) 
=- Lim 2212) =O. 
7->00 


Thus, u(z) iS a strictly convex decreasing function 
approaching zero as 2z27++t*? and hence u(z) > 0 for all gz. 


A sketch of u(z) 1S given in Fig. Aq-l. 


U 
O Z 
Fig. A-1l 
Smee uiz > 0 for all z, nz) = ou(z) > 0 for all z 
and therefore the diagonal term a (A + tn(z)) > 0 for 
QO 


melee Since A > 0 
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viz) = 54 6(2z) + amg? (z) - 2026(z)o(z) - 1(z) 
Define T(z) = 2¢°(z) - 2z9(z)¢(z) - 0%(z) , then 
V(z) = oe az ee art (Z) 


T(E =Azo- (2) + 22°6(z) Oz) + 226° (2) 
= 22° (z) (2) ~ eae 


= 22012) (290(zZ) — ¢(2)) 


= ~2z7¢6(z)u(z) (A.1) 
T(z) = -26(z)u(z) + 227$(z)ul(z) + 226(z) o(z) (A.2) 
Thus, 
Vi(z) =- = ZAguez, Jat Ta CZ) 
=- 2h OZ ete o (Zatz) L£rom (A.1) 
= -2296(2) ( = + TZ, ))) 


Gemee u(z) > 0 for all z 


> 0 LOLs Za<ae 
Vv  (Z) = 0 fOr. =o: 0 
< 0 £Or. Sze 60 


1 





Eee eS 2-4 (z) + TT" (2) 


V"(z) = - = = 
At z= 0 ’ 
T = 2A " 
v" (0) = - ace o(0) + wr" (0) 
= fa 6(0) — 276(0)u(0) from (A.2) 


Thus v"(0) < 0 and hence V(z) is maximized at z=0O and 


v(o) =  9(0) + mP(0) = 2A (.3989) + 7(.0683) >0 . 


lim V(z) = 0+7 lim T(z) =0 

Z>+0 Zeta 

lim V(z) = 0 +m lim T(z) = wW(-1) =-rT. 
Vg edo) T= Oo 


Peeretchn Of V(z) is given in Fig. A-2. 


V 





SS aLef EES & 


Since V(z) attains a maximum with positive value at z= 0 ana 
is decreasing for z > 0 and lim V(z) =0 , V(z) > 0 for all 


VAs th 
Ze 0, 
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